Dependence of the fluctuation-dissipation temperature on the choice of observable 
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On general grounds, a nonequilibrium temperature can be consistently defined from generalized 
fluctuation-dissipation relations only if it is independent of the observable considered. We argue that 
the dependence on the choice of observable generically occurs when the phase-space probability 
distribution is non-uniform on constant energy shells. We relate quantitatively this observable 
dependence to a fundamental characteristics of nonequilibrium systems, namely the Shannon entropy 
difference with respect to the equilibrium state with the same energy. This relation is illustrated on 
a mean-field model in contact with two heat baths at different temperatures. 
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Characterizing nonequilibrium states through general- 
ized, or effective, thermodynamic parameters is one of 
the important open issues in nonequilibrium statistical 
physics [H. One possible approach is to introduce ther- 
modynamic parameters conjugated to conserved quan- 
tities 0, H]. An alternative approach, more suitable 
for the definition of temperatures, is to generalize the 
fluctuation-dissipation relations (FDR) that relate the 
linear response to an external perturbing field with the 
correlation of spontaneous fluctuations [J]. At equilib- 
rium, the proportionality coefficient is precisely the tem- 
perature. In nonequilibrium situations, one may use the 
ratio between correlation and response as a definition of a 
non-equilibrium temperature, as proposed in the context 
of hydrodynamic turbulence , spin-glasses Q-lll , gran- 
ular materials [ol-flH. or sheared fluids [12] ■ The FDR 
approach is particularly useful since it allows for exper- 
imental 13 -3| and numerical (T^-HSl measurements. A 
necessary condition for a consistent definition of a FDR- 
based temperature is that its value does not depend on 
the choice of observables underlying the FDR. For glassy 
systems, it has been shown that no dependence on the ob- 
servable appears at the mean-field level 0, 0] . This prop- 
erty was also reported in numerical tests on more realistic 
models , althorigh the conclusions may depend on the 
model considered [8| . For non-glassy systems driven into 
a nonequilibrium steady-state, the situation remains un- 
clear, and no generic conclusion has been reached, even 
though a lot of work has been devoted to the study of 
generalized FDR [2^-[2^. Theoretical arguments on the 
observable dependence of the FDR-temperature are thus 
highly desirable. 

In this Letter, we study how such a dependence on 
the observable emerges in a specific class of stationary 
nonequilibrium systems. We study the time-dependent 
linear response of a family of observables to an exter- 
nal field. Relating these response functions to the as- 
sociated correlation functions provides us with a set of 
FDR and with the corresponding effective temperatures. 
These temperatures are found to depend on the observ- 



able, a property that we trace back to the non-uniformity 
of the phase-space distribution, measured with the Shan- 
non entropy. A quantitative relation between observable 
dependence and Shannon entropy difference with a refer- 
ence equilibrium state is obtained in a low forcing limit. 
We illustrate these results on a fully-connected model in 
contact with two heat baths at different temperatures. 

Considering a generic system described by a set of 
N variables Xi, i = 1, . . . ,iV, we introduce a family of 
observables Bp = J2iLi ^^'^^ ^ with p > an integer 
number. A small external field h, conjugated to an ob- 
servable M, can be applied to probe the system. The 
linear response of Bp to the external field is defined ac- 
cording to the following protocol. The field h takes a 
small non-zero value at times t < 0, and we assume that 
the steady state is established. At time t — the field 
is switched off and the time evolution of the observable 
Bp is recorded. The linear response Xp{i) is defined as 
Xp{t) = d{{Bp{t))) /dh\h=o, the average ((■■■)) being taken 
over the dynamics resulting from the above protocol. A 
FDR holds when the response Xp{t) is proportional to 
the correlation function (computed for h = 0) 



Cp{t) = {{Bp{t) - {Bp)) (Af(0) - (A/))) , (1) 



namely 



Xp{t) = ^ Cp{t) . (2) 



The proportionality factor is the inverse of the effective 
temperature Tp, which could a priori depend on p, and 
thus on the observable. Formally, {{Bp{t))) can be ex- 
pressed as 



N 



{{Bpim^ / l[dxMBpG'i{{x,}\{x'^)V{{x',},h) 

(3) 

with Bp = Bp{{xi}), and where Gf{{xi}\{x'^}), the zero- 
field propagator, denotes the conditional probability to 
be in a microstate {xi} = {xi,i = 1,...,N} at time t 



2 



given that the system was in a microstate {x[} at time 
t = 0, in the absence of the probe field. The distribution 
Vdx^}, h) is the stationary distribution of the microstate 
{x^} in the presence of the field h. Taking the derivative 
of Eq. ([3]) with respect to ft. at /i = 0, and using the 
relation dV/dh = VdhiV /dh, we get 



(4) 



the average being computed at zero field [30-1321] . 

To proceed further, an explicit form of the distribution 
Pdxilj/i) is required. As a general framework, we con- 
sider a class of stochastic markovian models, where a con- 
served energy E = X]t=i^ft.(^i) randomly exchanged 
between the internal degrees of freedom and with the 
environment (e. g. reservoirs or external forces). The ex- 
ternal sources and sinks drive the system into a nonequi- 
librium steady state. The resulting drive can be encom- 
passed by a dimensionless parameter 7 (e.g. a normal- 
ized temperature difference or external force). For zero 
driving, the dynamics satisfies detailed balance and the 
system is at equilibrium, characterized by the Gibbs dis- 
tribution Veq{{xi}, h) = Z^"" exp[-/3 Yl!i=\ ^h{xi)\ where 
/3 = 1/T is the inverse temperature imposed by an exter- 
nal bath, and Zj^ is the normalization factor. 

To simplify the calculations, we assume that the de- 
grees of freedom are statistically independent, namely 
V{\xi\, h) = Yl^^iP{xi, h). We thus focus on the single- 
variable probability distribution p(x, h). Considering the 
small driving limit I7I <ti 1, we expand the steady-state 
distribution p{x, h) around the equilibrium distribution 
Poq(a;, h) = ^ exp[-/3e,j(a;)] as 



p{x, h) = peq(x, h) [1 + iF{eh{x)) + 0(7^)] 



(5) 



The normalization of p{x,h) and Peq{x,h) imposes 



0, where 



is the equilibrium aver- 



age, and e denotes eh{x). Note that il p{x,h) follows 
Eq. ([S]), the factorized iV-body distribution V{{xi},h) 
is in general no longer a function of the total energy 
E = X^ili is thus not uniform over the 

shells of constant energy. As a result, the Shannon en- 
tropy of the nonequilibrium state should be lower than 
the entropy of the equilibrium state with the same en- 
ergy. Hence the entropy difference between the equilib- 
rium and nonequilibrium states with the same average 
energy provides an interesting characterization of the de- 
viation from equilibrium. Practically, the entropy differ- 
ence is determined as follows. We compute the average 
energy E{P, 7) of the out-of-equilibrium system, and we 
find the temperature (3* such that E{l3,-f) = Ecq{/3*), 
where E^qiP*) is the equilibrium energy at temperature 
(3*. From the factorization property of V{{xi},h), the 
Shannon entropy of the whole system is the sum of the 
entropies associated to each variables x^, so that we only 
need to compute the Shannon entropy per degree of free- 
dom S = — J dxp{x, h) lnp(x, h). The entropy difference 



AS* is then defined as AS = S'eq(/3*) - S'(/3,7), where 
Scq{f3*) is the equilibrium entropy at temperature /?*, 
and 5(/3,7) is the Shannon entropy of the nonequilib- 
rium state in the presence of a forcing 7. After some 
algebra, one finds 



AS* 




(6) 



Note that the 7^-term in the expansion (O of p{x, h) 
needs to be taken into account in the calculation, but 
eventually cancels out. It can be checked that AS" > 0, 
although this property is not explicit in Eq. The 
equality AS* = is obtained for a linear F{e), as in this 
case, p{x, h) can be recast into an equilibrium form with 
an effective temperature -see Eq. ([S]). Turning to the 
non-linear case, we parameterize F{e) as F{e) = a + 
be + ?7/(e), where 77 characterizes the amplitude of the 
nonlinearity. The parameter a is fixed by the constraint 
(F(e)) = 0. We then find A5 = j^rj'^a, where cr is a 
constant depending on the detailed shape of the functions 
/(e) and £^(2;). For instance, in the case /(e) = e^ with 



ft = and eo(a;) = ^x^, one has AS" = /P ■ 

We now come back to the FDR. As a simplifying 
hypothesis, we assume that the dynamics is such that 
each event decorrelates the involved variables Xi from 
their previous values (sec below for an explicit example) . 
Hence all correlation functions arc proportional to the 
persistence probability $(i), that is the probability that 
no event invol ving a given variable Xi occurred between 
times and t [2J]. Expanding for small field ft the lo- 
cal energy eh{x) = eo(a::) — hi(j{x) + 0{h?), with il^{x) 
an odd function, we find Cp{t) = N {x{tfP+^^{x{fi))) = 
N{x^P^^i(j{x)) $(t), the last average being performed on 
the steady-state distribution. A similar calculation start- 
ing from Eq. ^ yields for the response function 

Xp{t) ^ N[l3{x'P+^ij(x))^Yx^P+^^{x)F'{eo))] $(t) (7) 

where eo stands for eo(a;). To first order in 7, the average 
in the second term can be replaced by the equilibrium 
average. Expressing $(t) as a function of Cp{t), we find 
that the FDR |2]) is obeyed, with Pp = given by 



(x2p+V(a;)F'(eo))c 



(x2p+iV'(x)) 



(8) 



oq 



As expected, Tp generically depends on p, that is on the 
observable. Yet, a linear F(e), namely F{£) = a + be, 
yields an effective temperature fip = /3 — 76 that does not 
depend on the observable. A dependence on the observ- 
able arises when i^(e) has a nonlinear contribution. Using 
the parameterization F{e) = a + fee + ijf{e), we get from 
Eq. ^ that /3p — /3p=o is proportional to 77. Hence the 
dependence of the fluctuation-dissipation temperature on 
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the observable is directly related to the amplitude of the 
non-linearity in F{e). As this amplitude 77 is also cap- 
tured by the entropy difference AS", it is interesting to 
relate quantitatively /3p — /3q to AS*. We obtain 



(9) 



where Kp > is a dimcnsionless constant, depending on 
p, on the functional forms of /(e) and eii{x), but not on 
7 and rj. As an example, one finds in the case f{e) 



and eh{x) 



hx + 0{h'^) that Kp = Ap/V^. 



Thus it turns out that the dependence of the fluctuation- 
dissipation temperature on the choice of observable is a 
direct measure of the deviation from equilibrium. 

To illustrate this result, we consider an energy trans- 
port model on a fully connected geometry, in contact 
with two heat baths at inverse temperatures /3i and /32- 
The contact is characterized by a coupling constant v. 
A random fraction fi of the local energy Si = eh{xi) = 
I (xi — h)"^ is transferred from a site i to an arbitrary site j 
with a probability rate if{^\£i). After the transfer, Xi and 
Xj are changed according to x[ — h — ± ^ [xi — K)^ — 2fi 
and Xj~h — ±\/ {xj — h)'^ -f- 2/i with equiprobablc and 
uncorrelated random signs. The heat baths are char- 
acterized by injection rates Jq(/z) = ve~^°^ [a = 1,2), 
while transfers from the site i to any of the heat baths fol- 
low the rate i^(/3(/^|ei). We choose (^(// 1 Ei) = {l—^/si}^^/'^ 
to ensure that the equilibrium distribution is recov- 
ered when /3i = /32 [33(]. A sketch of the dynamics is 
shown on Fig. [T] Due to the fully connected geome- 
try of the model, the different sites become statistically 
independent in the thermodynamic limit iV — 00, so 
that the single-site distribution p{x, h) provides an ex- 
act description in this limit. In addition, as the dy- 
namics involves only redistributions of the local energy 
Ei = ^(xi — /i)^, the model can be effectively described 
in terms of local energies. The stationary distribution 
p(x, h) is then completely determined by the distribu- 
tion P{e). This change of variables is expressed in the 
probability densities as p{x,h) = ^\x — h\P{^{x — h)^). 
The equilibrium distribution at temperature T = (3^^, 
Pcq{x,h) = -y//?/27rcxp[— /3(a; — ^)^/2] corresponds to 
Pcqi^) = \fWh^ exp(— /?£). The master equation for 



other 

sites 
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FIG. 1: (Color online) Scheme of the mean-field model. A 
single site contains an amount of energy e. It is in contact with 
two baths at different temperatures Ti = fi^^ and T2 = fi^^ 
and with the other sites. 



the time-dependent iV-site probability distribution can 
then be recast into a nonlinear evolution equation for 
the single-site energy distribution P(e, t): 

dP , 



dt 



d^i [ Ji (/^) + J2{^i) + <^in(M, t)\P{£ ~ t) 

dfj, [Jlifi) + J2(m) + (Pin{^J., t)]P{£, t) 

+{2iy + l) / dfiip{fj,\£ + fi)P{e + fi,t) 
Jo 

-{2u+i) r dii^w)p{£,t) (10) 



where ipin{fJ.,t) = de (f{n\e)P{e,t) accounts for the 
energy transfers coming from all the other sites. 

To determine the steady-state distribution P{s), we 
consider the limit of a small temperature difference /3i = 
/3(1 — A) and /?2 — /3(1 + A), with A <C 1, and expand 
the distribution in A. The linear term in A vanishes be- 
cause the two heat baths play a symmetric role. The 
leading correction should thus behave as A^, so that the 
distribution can be written in a form similar to Eq. (|5|. 
namely 



P(e)-Peq(£) [l + \^F{e)+0{X')\ 



(11) 



This scaling form is validated by direct numerical simu- 
lations of the stochastic dynamics, as shown in Fig. [21[a). 
Simulations were performed on a system of size N = 102, 
with (3 = 1. The relatively small system size allows for 
long time averaging until a time T^ax ~ 10'' or 10^, so 
as to reach a satisfactory statistics (in a unit of time, all 
sites have in average experienced about one redistribu- 
tion event). In Fig.[5Jb), the function F{e) obtained from 
simulations is shown for different values of the coupling 
strength v. We observe that the curvature of F{e) is re- 
duced when decreasing v. Eq. PH)) has no exact solution 
involving a finite polynomial function F{e). To flnd the 
best polynomial approximation F^^''{e) = 'Yl^k^a^kP^s^ 
at a given order i, we devised a variational procedure {j3 
factors are introduced to make Ofc dimcnsionless). The 
function _F'-^^(e) is obtained analytically by minimizing 
the error, under the constraints of normalization and zero 
neat flux with the baths, in the evolution equation (|10p 
linearized in F{e). The error is defined as the equilibrium 
average of the square of the r.h.s. in the linearized equa- 
tion. For L > 2, we find that the coefficients a^, fc > 2, 
in the expansion are numerically small, as illustrated in 
Fig. EJc) . A second order polynomial is thus already a 
good approximation of F{e) for < 1 [see Fig. [l{b)]. 
Taking into account higher order terms in i^(e), we find 
that the relation © between the observable dependence 
and the entropy difference is valid to a good accuracy 
[Fig. Efd)]. In addition, we observe that in the limit 
ly 0, the coefficients a^, k > 1 vanish while Oq — > — | 
and fli — >■ |. Therefore the temperature Tp becomes ob- 
servable independent in the small coupling limit. 
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FIG. 2: (Color online) (a) (P - Po)/A^ as a function of the 
local energy e obtained by numerical simulations of the dy- 
namics {u = 1, Tmax = 10®). (b) The function F(e) obtained 
by simulations (noisy lines) in comparison with the analyt- 
ically obtained results for F^^\s) (solid lines) for different 
values of (A = 0.2, Tmax = 10^). (c) //-dependence of the 
coefflcients in F^^\e) (see text), (d) Parametric plot in v 
of A/3p = \l3p - l3o\/l3 versus \/AS obtained using F'^'(e), 
either fitted to numerical data (x) or calculated in the ana- 
lytical approximation (solid lines). Dashed lines: Eq. ((9| with 
Kp — 4p/\/3. Simulation parameters: A — 0.05, Tmax = 10^. 

To sum up, we have shown that a large class of 
nonequilibrium systems generically exhibit observable 
dependence of the fluctuation-dissipation ratio, even in 
the mean- field case. Accordingly, a unique nonequilib- 
rium temperature cannot be defined from the FDR. The 
dependence on the observable can be traced back to the 
non-uniformity of the phase-space distribution on shells 
of constant energy, quantified by the difference of Shan- 
non entropy between the equilibrium and nonequilibrium 
states with the same energy. We have illustrated these 
results explicitly on a mean-field model connected to two 
heat baths, confirming that observable dependence ap- 
pears in the driven stationary state. The dependence 
however becomes weaker when the coupling to the reser- 
voirs is decreased. This might be the reason why the 
observable dependence has not been encountered in nu- 
merical simulations of granular gases [2^. Furthermore, 
it turns out that the entropy difference AS" is a relevant 
characterization of nonequilibrium systems. If AS' = 0, 
a single temperature emerges from the FDR, and the 
statistical properties resemble closely that of equilibrium 
systems. In contrast, if AS* > 0, the system can be de- 
scribed by two parameters, a reference temperature (e.g., 
Tp^o) and AS. It would be interesting to try to mea- 
sure AS experimentally or numerically through the use 
of the FDR, in real out-of-equilibrium systems like gran- 
ular gases 23 1 or turbulent flows (s^ . 
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